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Abstract 

This paper presents a systematic investigation of the integrability conditions for 
nonautonomous quad- graph maps, using the Lax pair approach, the ultra-local sin- 
gularity confinement criterion and direct construction of conservation laws. We show 
that the integrability conditions derived from each of the methods are the one and 
the same, suggesting that there exists a deep connection between these techniques 
for partial difference equations. 
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tions; conservation laws 



1 Introduction 



Recently, there has been rapid development of research into integrable discrete 
nonlinear systems governed by both ordinary difference equations (mappings) 
and partial difference equations (lattice equations) [T|2|3|Hl5] . It has been sug- 
gested that discrete systems governed by difference equations are more funda- 
mental than the continuous ones described by differential equations. Efforts 
have been made by several research groups to develop analytical techniques 
to determine whether or not a given nonlinear partial difference equation is 
integrable [6. 7. 811^ . It has been possible to identify several rich classes of in- 
tegrable nonlinear autonomous partial difference equations that are amenable 
to solution by linear methods [9|10|ll|12lll3lll4fl5|16j . For a given autonomous 



Partially supported by the Indian National Science Academy and the Royal So- 
ciety. 



Preprint submitted to Elsevier 



February 8, 2008 



integrable nonlinear partial difference equation, the question of finding an in- 
tegrable nonautonomous version systematically is not yet solved. Examples of 
nonautonomous quad-graph equations are considered in P^PTj : in particular 
[T7] deals with reductions to ordinary difference equations. 

The aim of this article is to demonstrate that the combination of the Lax pair 
formalism, the ultra-local singularity confinement criterion and construction 
of conservation laws provides an efficient tool to investigate integrable nonau- 
tonomous nonlinear partial difference equations. We use these methods to find 
integrable nonautonomous versions of the discrete Korteweg-de Vries (dKdV) 
and modified discrete Korteweg-de Vries (dmKdV) equations. 

The plan of the article is as follows. In §2 we explain the Lax pair formalism 
for nonlinear partial difference equations on the quad-graph and derive con- 
ditions for the existence of a Lax pair for nonautonomous dKdV and dmKdV 
equations. In §3 we derive conditions for each of the above nonautonomous lat- 
tice equations under which the ultra-local singularity confinement criterion is 
satisfied. Section 4 deals with the derivation of conditions on the above nonau- 
tonomous lattice equations for which three-point conservation laws exist. In 
§5 we draw conclusions and describe a possible direction of research. 



2 Lax pair formalism for partial difference equations 



In the following, uf^ denotes the value of the dependent variable u at the point 
(/,m) e Z^. Attention is restricted to difference equations on the quad-graph, 
which are equations of the form 

We also assume that this expression can be solved for any one of the values 
of u; in particular, we write 

uJlf = u;{l,m,uT,uT^„ur'). (1) 

Such a nonlinear partial difference equation is said to be completely integrable 
if it arises from the compatibility condition of a system of linear difference 
equations: 





= L{1, m; k) 


vrik) 




'vr\k)' 






wr\k)_ 
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M{l,m; k) 



vrik) 

wr{k) 



(2) 

Here f[" and are wave functions defined on the quad-graph (or at the 
nodes of a two-dimensional lattice) as functions of a spectral parameter k. 
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The 2x2 Lax matrices L{l,m; k) and M{l,m; k) describe the change in the 
wave functions under a horizontal and vertical shift respectively. They depend 
upon the spectral parameter k and on u, which plays the role of a potential. 
As ([2]) gives two different ways to express v^^{k) and w]^^{k) in terms of 
v["(/c) and wY^{k), this leads to the compatibility condition 

M{1 + 1, m; k)L{l, m; k) - L{1, m + 1; k)M{l, m; k) = 0. (3) 

We assume that L{1, m; k) and M(/, m; k) depend on the potential only through 
{u^.n^i) and [uY^^u^^^] respectively, and that the compatibility condition 
([3]) yields the integrable quad-graph equation ([1]). 



2.1 Lax pair compatibility conditions for the nonautonomous dKdV equation 



Generalizing the Lax pair for the dKdV equation, we consider Lax matrices 
L(/, m; k) and M(/, m; k) of the form 



L(/, m; k) - 
M{l,m; k) 



u 



1 



m+l 



U 



m+l 
I 



e + (3r + - uf+'){dT + uT) dr + < 



m+l 



(4) 



At present, we regard a, f3, a, b, c, d as arbitrary functions of / and m that may 
occur in the nonautonomous dKdV equation. The compatibility condition ([3]) 
for the Lax matrices (jl]) has four components that constrain the arbitrary 
functions. 



Components (1,1) and (1,2) of the compatibility condition give 



u 



m+l 
l+l 



l,m+l 
"Z+1 



a. 



u 



m+l 
/ 



U 



l+l 



Oi -t- a^+i 



(5) 



and 



aT^f + dT = cT + bT^f. (6) 

Substituting these results into the (2,1) component of the compatibility con- 
dition yields 



A+i = PT. 



a 



m+l 



a 



I 1 



(7) 



^m+l I ^m+1 



jm I m 
"«+l ^ "i+1- 



Finally the (2,2) component produces no additional constraints. Equations ([7]) 
show that 

«z"^ = /(0, f3r = 9{m), (9) 
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for arbitrary functions / and g. To proceed further we need to solve ([6]) and 
(IH]). To do this we define new variables 0™ and ip]^ by 

0r = «r-fer, cr = cr-c^r, (lo) 
V'z+i = &r+cr, i^r' = dT+ar. (n) 

Then ([6]) and ([8]) imply that 
Let A™ be a solution of 

\m \m+l im 

^/ - ^i-i — 9i ) 

then 

- v^r^i = Ar - Art' - Ai^t^ + Ar'. (12) 

Note that A™ is defined up to an arbitrary function of / + m. The solution of 
(Ha) is 

^r = Ar'-Ari+M^+H, 

where h{l+m) is an arbitrary function. We can eliminate h{l + m) by replacing 
AJ" by A™ + H{1 + m), where H{1 + m) is a solution of 

H{1 + m + 1) - H{1 + m - 1) + h{l + m) = 0. 

Finally ([5]) can be rewritten as 

{uT^f -uT-bT^\'+bT- \T+\' + XT) («r ' - <i - &r ' + + A^i - Ar ^ ) 

= f{l)-g{m). (13) 

The transformation 

mJ" + 6;" + a;" 

maps f|T3|) into the standard dKdV equation 

(«-t' - OKU - O') = /(O - 9im). (14) 

Thus we find that all nonautonomous dKdV equations which admit the Lax 
pair representation (jl]) can be mapped into ([T^ by a point transformation. 
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2.2 Lax 'pair compatibility conditions for the nonautonomous dniKdV equa- 
tion 



By analogy with the autonomous dmKdV equation, let 



L{1, m; k) 



M{l,m; k) 



k^cV 



u 



m „ m 
I "'l+l 
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U 



m+1 



U 



m+1 



U, 



(15) 



Here a, b, c, d,p, q, r and s are arbitrary functions of / and m that may occur 
in the dmKdV equation. The (2,1) and (1,2) components of the compatibility 
condition ([3]) give 
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(16) 
(17) 



H+l\ 



Then the (1,1) and (2,2) components of the compatibility condition yield the 
constraints 

„m ^m+1 m ^m+1 

ft+1 _ Pi ^l+l _ 'l ^-j^g^ 



„m „m ' „m ».m 

HI Pi ^l 'l 

rjm+1 m jm m 

^ _ ^ ";+i _ ^ 

„m ^m ' m+1 hm 

H+1 ^l 



(19) 



Combining (|T6l) . ( ITTI) and ( ITSi) gives one further constraint on the unknown 
functions: 

um „m l,m+l „m+l 

„m ^m T-^+l ' ^ ' 

P« 'l Pi ^l 

In Appendix A, we solve the constraints f|T8l) . f|T9|) and fl20|) . obtaining a 
nonautonomous version of the dmKdV equation. This can be mapped by a 
point transformation into the standard dmKdV equation: 

m+i_^m[9im)ur^i-fii)ur'] 



[9{m)ur' - /(/)<+i] 
Here f{l) and 5'(m) are arbitrary functions. 
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3 The ultra-local singularity confinement criterion 



A difference equation is said to passes tlie singularity confinement criterion 
(which can be viewed as a discrete version of the Painleve Property) if every 
singularity is confined to a finite number of iterations [5]|llfl4j . Suppose that, 
in solving an initial- value problem for ([1]), a singularity occurs at the point 
(/, m), that is 

uY' = 0{e) or ur = oQ^ as e 0. (21) 
The ultra-local singularity confinement criterion requires that 



= 0(1), = 0(1), 

uT_f = o{i), «ri' = o(i). 



We can investigate the consequences of this criterion by solving the quad- 
graph equation for each of the variables in turn and applying the shift 
operators appropriately. 



3.1 Ultra-local singularity confinement conditions for a nonautonomous dKdV 
equation 



Consider a nonautonomous version of the dKdV equation that is of the form 

{uT^,' -uT+ BT) «i - +Ar) = cr, m) 

where A™, Bf^ and C[" are arbitrary functions. This is consistent with the 
generalization ([5]) that arises from the Lax pair. Singularity confinement with 
= 0(e) gives no constraints on these functions. However, if mJ" = 0(l/e), 
nontrivial conditions are obtained, as follows. Equation fl22|) can be written as 

= uT - BY" + ' 



H+l -^Z „,m+l I Am ' 

^Z+1 "~ ^Z + ^Z 
/^m—l 

^Z ^Z-1 + ^Z-1 



C 

m-l _ m _ 4^-1 _|_ I 
^Z+1 — ^Z ^Z + 



m— 1 



ut\' = uT + At, - 



rim 

'-'Z-1 
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Therefore 



ut\' = 0{l) ur'-nt^ = -At\'+0{e), 
uiri' = Oil) uYU-uT-' = -Br'+0{e), 
utV = 0(1) - uZ, = -BYU + 0(e). 



From these conditions we obtain 

At\'-AT = B\l,-BT-\ 
Equation (122|) can also be rewritten as 



m+1 _„,m _ nm j 



/-mi— I 

or 



m+1 _„,m _|_ 4' 



u^:,V -uT + BY 

Combining with ^ gives us 



+ 



or 



m+1 _ m _|_ 
^Z+1 ^/ + 



If = 0(l/e) then this yields 

w| = -(o^i + or' - on^ - cn + o(i), 

e 

where is any of uT_\\ u\Xi\ uTS^}. Therefore 

_j_ ^m— 1 ^m— 1 (J^ Q 

The general solutions of ([23]) and (12^ are 

Am \m \m+l 

A — ^« + l ~ ^« 5 

Oz™ = /(/)-^7(m), 
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where AJ", h{l + m), /(/), g{m) are arbitrary functions. After substitution of 
(IM]) into ([22]) we obtain 

_ uT + x^i\' - \T + h{i + m))(t.™i - + Aril - Ar^) 

= /(0-^M- (33) 

The transformation 

-XT + H{l + m), 

where H{1 + m) is the solution of 

H{l + m + 2) - H{1 + m) + h{l + m) = 0, 

maps ( l33l) into the standard dKdV 

- - ur') = f{l) - g{m). (34) 

Thus every nonautonomous dKdV equation of the form fl22l) that satisfies the 
ultra-local singularity confinement criterion can be transformed into (!34|) . 



3.2 Ultra-local singularity confinement conditions for the dmKdV equation 



We consider the following nonautonomous generalization of the dmKdV equa- 
tion: 



m+1 

where AT, BT and C™ are arbitrary functions. This equation can be written 



as 



m+l _ m ] _ Am m _ r>m„,m+l 

Ul+l [Ui 1^1 M^+iJ — Ui [/ii Ui^^ t^i 



U 



ut\'[At\'uT-' - BTrM-i] = <[uTi - cr>T-'] 
uT^~,'[AT-'uT-' + cr'uT^i] = uTK^, + Br'u^-'' 



,,m+lr m+l I r>m m ] _ m\ Am m , „,m+ll 

If = 0(e) then 

«^ = 0(1) ^ ;^ = ^ + 0(e), 

m-l R™-! 

<T^ = 0(1) ^ ;^ = -^ + 0(6), 
<V = 0(1) ^ ^ = -5-1 + 0(e). 
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Consequently 



If = 0(l/e) then 



0. (36) 



u, 



JDr 



m 



u 



uf_f = 0(l) ^ '^ = -tlt± + Oit). 



Therefore 



^-^^^ = 0. (37) 



AT Br cr 



Thus the nonautonomous dniKdV meets the ultra-local singularity confine- 
ment criterion only if fl36|) and fl37j) are satisfied. By solving these equations 
we obtain a nonautonomous version of the dmKdV equation that satisfies the 
ultra-local singularity confinement conditions (for details see Appendix B). 
This equation can be transformed to the standard dmKdV equation, 



where /(Z) and g{m) are arbitrary functions. 



4 Conservation laws 



In this section, we find conservation laws for quad-graph equations by the 
direct method, which is explained fully in p^|IT6] . Three-point conservation 
laws have components 

that satisfy the following functional equation on solutions of the given quad- 
graph equation: 

Fil + 1, m, „ <|i) - F(/, m, ^) + G(/, m + 1, <|^) 

-G{l,m,uT,ur)=^- (38) 
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To find these conservation laws, we first substitute 



into (138|1 . The resulting equation involves m™, u^li and m™^^, but each instance 
of F and G depends on only two continuous arguments. Therefore we can 
eliminate terms by repeated differentiation until a partial differential equation 
for F is obtained. Having solved that, it is simple to work up the hierarchy of 
functional-differential consequences of fl38l) until the general solution has been 
obtained. The same process can also be used to find higher conservation laws. 
However, the existence of three-point conservation laws is sufficient to classify 
the nonautonomous dKdV and dmKdV equations. 

The calculations are immensely complicated, and computer algebra is essen- 
tial; we have used the computer algebra system MAPLE. For brevity, we sketch 
the results of the computations, without giving full details. 



4.1 Conservation laws for the nonautonomous dKdV 



The generalization of the nonautonomous dKdV equation fl22p can be trans- 
formed to 

V^i+i + Bi ){Ui -Ui^^)-Ci , {69} 

by the point transformation 

where Q]" is a solution of 

QT+i - QT"-' + = 0, 

and where 

Br = Br + QT^i' - QT- 

This simplification greatly speeds up the computations, without affecting the 
number of independent conservation laws that exist. When the direct method 
is used to reduce fl5I?l) to a partial differential equation, we find that 

F =8 {urfuf+^cr 

+4nj"+^ {ur^^{urfvr +2ur+\r^r +'^f^rur - ^cr^^cr +2cr^^c,r^^) 
"+ir (i/™c7™+i - vrcr + 2 cr^^nf^+^ + 2 - 2 /x™) ur{ur+^fT'^ 



u 



I 2 \tt t L I I i It/ 



-4ni!ti (2cr'Q"^ + + uT+i^rwT? -^^^rcrur +2ur+,cr^r) 



- 2- 



cr (CD' 
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where T = vJ^C^ - i^r^r^^ - 2 /xf^ - 2 Here ^P./^F. ^F" and are 

functions which satisfy the following constraints: 

^ _ 4/i, _ (42) 

Note: at this stage, we have not completed the direct method calculation of 
the conservation laws, but the above necessary conditions lead to a substantial 
further simplification of the problem. 

The general solution of the system (HOi) is 

Cr = ^(/ + m), i^ = -E(l^m-\), z/r = ci(-l)'+™, (43) 

where H{l + m) is an arbitrary nonzero function and Ci is an arbitrary nonzero 
constant. Combining these results with fHTl) . and fH2|) . we obtain 

= -(-iy+''UH(l + m + l) - H(l + m- 1)). 

Cl 

Therefore three-point conservation laws exist only if the nonautonomous dKdV 
equation is of the form 

(uT^f -uT + |(-l)'+™[/J(/ + m + l)-H{l + m- 1)]) («r'-«l+i) = CT- 

(44) 

This equation is mapped by the point transformation 

uT^uT- -(-!)'+"//(/ + m - 1) 

Cl 

to 

«t^-«n(«r'-«m) = cr- (45) 

Therefore it is enough to seek conservation laws of ( l45l) . Applying the full 
direct method to ( l45l) gives us one further condition on C™: 

^m+l ^m+1 _|_ g 

Consequently all nonautonomous dKdV equations that have nontrivial con- 
servation laws can be mapped to 

«r - uT){u1l, - ur') = /(/) - g{m), (46) 
whose three-point conservation laws are are 



;i) F = (-1)'+"^ (^2<m5"+^ + g{m 

G = _(-!)'+- {2urur+, + /(/)) 
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(2) 
(3) 
(4) 



F 
G 
F 
G 
F 
G 



- \u 



u 



m+l\ 



\l+m 



m+1 
I 



{g{m))' 



4-2 Conservation laws for the nonautonomous dmKdV equation 



In the same way (but with fewer details), we apply the direct method to the 
nonautonomous dmKdV equation (1351) . Then the components F and G are of 
the form 

Am „ m T3m „ m . m 

Q / vJ^Ut"^ + _|_ j^'>n I I r^l 

' TDm I '+1 ' „,m Am „.m /^m„,m„,m ' 

where the functions ^, fi, v and C satisfy the constraints 



tm /-m 


— 'il ' 


m m 
H'l+l'^l+l 


= f^>r. 




tm / tm+l 
— Si+l/ Si 


AT 


^m / fm 
~ Si+l/ Si ; 


cm+1 m+1 


m /-m 

— H'l ! 


m+1 Am+l 


— '^l ' 


,,m Tjm 


m/^m 
— '^l ■ 



The general solution of these constraints is very messy. However, it yields the 
result that the only nonautonomous dmKdV equations with nonzero A, B, C 
that admit conservation laws can be transformed to the standard dmKdV 
equation 



m+1 ^ ^m 



The three-point conservation laws for this equation are 

m m+1 

(1) F = ^^i^, 
g{m) 

f{l) ' 
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5 Conclusion 



We have studied tlie existence of nonautonomous versions of dKdV and dmKdV 
that satisfy the Lax pair constraints, the ultra-local singularity confinement 
criterion and the conditions under which a quad-graph admits three-point 
conservation laws. Each of these conditions imply that all integrable systems 
of these forms are related to the standard dKdV equation 

- uT) «i - ur') = f{l) - g{m), 



or dmKdV equation 



respectively, by a transformation 

uT ^ P{l,m)uT + Q{l,m), 

where P {1,171) and Q{l,m) are arbitrary functions. This result is in agreement 
with results from the paper [9] which uses the consistency approach. It suggests 
that there is a deep connection between the various criteria for integrability 
of quad-graphs. In order to find this connection more research has to be done. 

The methods that we have used in this paper can be applied to the other 
integrable quad-graph equations in the classification by Adler, Bobenko & 
Suris 0. These equations all have the tetrahedron property; at present, it 
is not yet known whether there exist integrable quad-graphs without this 
property that are not linearizable (see [T8lfl9] ). If such quad-graphs exist, 
they could also be tested by the methods in this paper. 
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APPENDICES 



A Lax pair compatibility conditions for the nonautonomous dmKdV 
equation 



By combining equations (fT9l) . we obtain 



l+l 



This can be integrated (assuming that the domain has trivial difference coho- 
mology) to give 

(A.l) 

where ?/^[" is an arbitrary function. Similarly from f|T8l) we have 

m+l rA'^+l 

where 05" and xT arbitrary functions. Equation ( l20l) also can be integrated: 

^ = i^(0, (A.5) 

'l 

where H{1) is an arbitrary function. By combining ( lA.ll) . ( 1A.3I) and (lA.Sp we 
obtain 



(A.6) 

Let H{1) = where is defined up to an arbitrary nonzero constant 

factor. Then 



C..Xiy(' + l) ,gXf^_ and.oC = Cxrf(OGM. (A.7) 
Therefore, from (lA.ip . 



(A.8) 
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and, from (lA.2p . 



m+l 



r'xr'Gjm + 1) 

cr^rxrG{m) 



(A.9) 



Substituting ( lA.SI) and ( ]A.9I) into the first of the remaining constraints ( |T9 
gives 



= [hT^'<Pr\r'F{l)G{m + l)]\hT(t^TxTF{l)G{m)]. 



Solving this equation, we obtain 



m m+l 

'h '11 



^■•pxTF{l)G{m) ' 



d 



' 4>TxTF{i)G{m) ' 

<Pr\r'F{l)G{m + 1) 



where ?7™ is an arbitrary function. Therefore flTB]) amounts to 



where {f{l)f 



T.+1 _ Vl+1 Xl '^1 



TT 



m + l -v/^ 



v 



m+l ""^ 



M, - 



F(0 



^ and (5'(m))^ 



G(m) 
G(m+1) 



. The transformation 



(A.IO) 
(A.ll) 



(A.12) 



where ^'^"^^V 
equation 



Xl 



g{m) and -4^^ = /(O; reduces flA.12p to the standard dmKdV 



u 



m+l 

i+i 



u 



\g{m)uT^,-f{l)ur'] 



B Ultra-local singularity confinement conditions for the nonau- 
tonomous dmKdV equation 



The conditions fl36|) and fl37j) amount to 



Consequently 



Am Tjm—l/^m jDm/^m+l ' 0^'^) 



Tym—l/^m 
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where H{1) is an arbitrary function. Integrating this gives 

5- = (/(/))2!^ Cr = ^, where f/iii-llV = 
Therefore f IB.ip yields 

m— 1„ 



and so 



where g{m) is an arbitrary nonzero function. Then ([35D amounts to 



'+1 „/,m I „l,rn„,m+l „/,ni-l„,m 



The transformation 



where = 5'(m) and -^^(jy^ = /(O; reduces ( 135!) to the standard dmKdV 

equation 



\g{m)uT^,-f{l)ur'] 
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